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Abstract 

A comparison between an analytical calculation of the polarizability of a 
mesoscopic interacting system in the random phase approximation and numer- 
ical exact diagonalization results is presented. While for weak interactions the 
analytical calculation fits the numerical results rather well, deviations appear 
for stronger interactions. This is the result of the appearance of intermediate 
range correlations in the electron density, which suppresses the polarizability 

below its classical value. The relevance to quantum dot systems is discussed. 
PACS numbers: 71.55.Jv,73.20.Dx,71.27.+a 
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The study of the polarizability of small mesoscopic systems has a long history. Already 
in 1965 Gorkov and Eliashberg ]I] have attempted to calculate the average polarizability of 
a small metallic grain by applying the concept of energy level repulsion from the random 
matrix theory (RMT) Their calculation, which can be backed by a non-linear sigma 

model calculation [[|], shows that the average polarizability of a mesoscopic system (a) oc 
(nL)^ 1 ^ L 3 (here L is the length of the system and the inverse screening length k^" 1 ) = 
Sde 2 N(0), where N(0) is the level density at the Fermi energy, Sd = 2ir, (4n) for two (three) 
dimensional systems and (...) denotes an ensemble average). This result is much larger than 
the classical value of the polarizability a a L 3 . The reason for this apparent enhancement 
of the polarizability is, as pointed out by Strassler et. al. ||, the neglect of electron-electron 
(e-e) interactions which reduce the quantum corrections dramatically. Once kL ^> 1, which 
for metallic systems corresponds to L ^> 1A, the polarizability is proportional to the volume 
||. The effect of e-e interactions in the random phase approximation (RPA) has been 
recently incorporated into the RMT formulation of the polarizability by Efetov [J]. The 
fluctuations in the grand-canonical ensemble were shown to be strongly suppressed by the 
e-e interactions , while the suppression is less significant in the canonical ensemble []TU 



In recent experiments on disordered quantum-dots [O it turned out that the ground 



state energies of those dots show large fluctuations |[l^J13[| . These fluctuations are larger 
than expected from RMT theory in which the e-e interactions are incorporated by the RPA 
13| . This deviation from RMT theory is the result of the appearance of intermediate range 



correlations in the electron density as a result of strong e-e interactions. In this paper we 
shall show that the average polarizability is also affected by strong e-e interactions. While 
in the weak interaction regime the average polarizability is proportional to L 3 , for strong 
interactions the polarizability tends to zero. 
The polarizability may be written as 

a = e 2 ^— ( df(z - L/2)n(r)$(r), (1) 
oE J 

where n(r) is the electronic density, $(r) is the local electrostatic potential as a result of 



applying an external electric field Ez and the polarizability is calculated in the z direction. 
We shall begin by considering the weak interaction limit. In this limit one can treat both 
the disorder and the e-e interactions perturbatively. Taking into account the e-e interactions 
in the RPA one may rewrite Eq. ([I]) for the average polarizability as || 

(a) = e 2 N(0) J dfdff(z - L/2)x{r,r/)(zf - L/2), (2) 

which after assuming a rectangular geometry and performing a Fourier transform results in 

/ T\2d 

(«) = e 2 iV(0)(-j J2[z(q)] 2 X (q), (3) 

where q = (n x ir/L)x + (n y ir/L)y + (n z 7f/L)z with n x ,n y ,n z = 0,1,2... (of course for 2D 
systems q = (n x n/L)x + (n z n/L)z), 

z(q) = T 5 qx=0 5 qy=0 — , (4) 

-k q z 

and 

X{ql = (l) 1 + N(Q)VU > (5) 

for short range interactions represented by an interaction potential U(f,ff) = VUS(f — fl) 

(V = a d where a is the range of the interaction and d is the systems dimensionality). For 

the Coulomb interaction U(r,rf) = Ua/\f — ff\, where U is the strength of the Coulomb 

interaction between two particles at distance a, and 

[2\ d q( d ~ l ) 
X ^ = [lJ q( d -V + S d N(0)alI ' (6) 



Inserting all the definitions in Eqs. (|4]-|6|) and performing the summation in Eq. (Q) results, 
for short range interactions, in 

where N(0) = (L A)" 1 , and A is the single electron level spacing. For long range interactions 



Note that by using the constant U to describe the strength of interactions between the elec- 
trons while using e to describe the interaction between the electric field and the electrons 
we artificially create two different interaction scales. Nevertheless, this distinction is useful 
because it enables us to treat the interactions between the electrons by some appropriate 
effective interaction (say, short range interactions) while the interaction between the elec- 
tric field and the electrons contributes a prefactor (kL)^ 1 ^ to the polarizability. Only in 
this sense can one understand the original treatment of Gorkov and Eliashberg |l| of the 
polarizability of non-interacting electrons. If we insert U = e 2 /a in Eq. @ we obtain 

<«> = £ + {kLMW))' 1 , (9) 

which for kL 3> 1 gives the expected classical behavior 

~ ^r L >- < 10 > 

Thus the RPA approximation corresponds to the classical result a cc L 3 . 

In order to confirm the above results and to check their range of validity we have per- 
formed a numerical calculation of the polarizability for a system of interacting electrons 
on a 2D cylinder of circumference L x and height L z . We used the following tight-binding 
Hamiltonian: 

H = £( e fcj + (j ~ {L z + l )/ 2 ) E )al d a kJ - Vj2( a l,j+i a k,j + 4+i,i°fcJ + M + H in t (H) 

k,j k,j 

where a k is the fermionic creation operator, e kt j is the energy of a site (k,j), which is chosen 
randomly between —W/2 and W/2 with a uniform probability and V is a constant hopping 
matrix element. We set the coupling the electrons and the external electric field E to be 
unity (i.e., e = 1). H int is the interaction part of the Hamiltonian which for the short range 
interactions is given by: 

Hint = U a kJ a kd a l,p a l,Pf ( 12 ) 

{k,j>l,p} 

where {. . .} denotes the nearest-neighbor pairs of sites, and for the Coulomb interaction is 
equal to 



Hml = v y, (is) 

k,j>i, P l rfc .i r i,p\/ s 

where K is a positive background maintaining the charge neutrality and s is the lattice 
constant. 

For a sample of M sites and N electrons, the number of eigenvectors spanning the many 
body Hilbert space is m = The many-body Hamiltonian may be represented by an 

m x m matrix which is numerically diagonalized and the ground state eigenvector ty E (k,j) 
is obtained for different values of the e-e interaction and electric field E. Here we consider 
a 4 x 4 lattice (i.e., M — 16 sites) and N = 4, 8 electrons. We chose W = 8V for which this 
system is in the metallic regime |l4j and average the results over 500 realizations for each 
value of interaction strength. The polarizability is calculated by 

a=m -w) 

E 

where 

d(E)=J2\^ E (k,j)\ 2 (j-(L z + l)/2). (15) 

k,j 

Our main goal is to study whether the effect of the e-e interactions is predicted correctly 
by Eqs. (0-(^). Therefore it is convenient to plot (a(U = 0)) / (a(U)) as function of interaction 
strength, which for short range interactions are is shown in Fig. [l]. In order to evaluate Eq. 
(0) for the lattice case one should replace V/L d by Zj [M — 1) where the average number of 
nearest neighbors, Z, in our case is given by Z = 3.5 due to the presence of boundaries. The 
average level spacing is A = 0.58V resulting in (a(U = 0)) / (a(U)) = 1 + 0AU plotted in 
the figure. We can see a reasonable correspondence to the numerical values for both values 
of N for weak interactions, i.e., U < V. For stronger values of interactions strong deviations 
appear for N = 8, while for N = 4 the correspondence holds. We shall discuss the reason 
for this behavior shortly. 

For the long range interactions the results are shown in Fig. For a lattice system 
one should replace S 2 = L^ 1 J d 2 r/r = 2tt in Eq. (|8|) by S 2 * = L^ 1 J2k,j^i, P \^k,j — 



which in our case results in S2* = 2.44. Thus, after incorporating S2* in Eq. @ one obtains 
(a(U = 0))/ (a(U)) = 1 + 0.32U plotted in the figure. Again a good agreement between the 
theory and the numerical results is seen for weak interactions (U < 4V), while deviations 
appear for both values of N for strong interactions. 

Thus, the RPA describes the polarizability rather well as long as the interactions are not 
too strong, resulting in the classical polarizability of Eq. (|l0j) a oc L 3 . Once the interactions 
are strong the polarizability is suppressed compared to the classical value, i.e., a<L 3 . 

The reason for the deviations is the appearance of intermediate range correlations for 
strong interactions. This may be clearly seen by defining a two point correlation function: 

_ ^k,j>i, p C(r ktj — ri^)5\ff k ^_ff hp \ %r ^ 

J2k,j>l,p d\r kJ -r^ p \,r 

where 

c{ k ' 3 " l ' p) ~ \ m^mmm / • (17) 

In Fig. |3| we present the correlation C(r = y/Es) for different values of interaction. Under 
the RMT assumptions, as well as for RPA, we expect C(r = y/Es) — > 0. On the other 
hand, once Wigner crystallization occurs we expect C(r = y/Es) — > — 1 for N = 8, while 
C(r = Vbs) -> 1 for iV = 4. It can be clearly seen that for weak interactions (i.e., U < V 
for short range interactions and U < AV for Coulomb interactions) C(r — y/Es) ~ 0. For 
stronger interactions a transition towards the Wigner crystal values of C(r = Vbs) is clearly 
seen for both the Coulomb and the short range interactions in the half filled case (N = 8). 
This correlation is a sign of the appearance of intermediate range order in the electronic 
system. We call this intermediate range order and not long range order since it decreases 
on the order of several lattice sites. Only for stronger interactions does a real, long ranged, 
Wigner crystal appear. Naturally, for lower fillings (N = 4) the transition is much weaker 
for the Coulomb interactions and absent for short range interactions. 

This behavior is the explanation for the deviations from the perturbative results of the 
polarizability at strong interactions seen in Figs. As long as there is no interme- 

diate range order, RPA holds. Once intermediate range order appears, the electrons are 



correlated, i.e., the electrons have a spatial order and exhibit a strong resistance to be- 
ing compressed. Therefore, the application of an external electric field has less influence 
on their spatial distribution and (a(U » V)) — > resulting in the strong deviation of 
(a(U = 0))/(a(U » V)) from the linear dependence especially for the half-filled cases in 
which the effect of intermediate range order is the largest. For N = 4 the deviations towards 
larger values of (a(U = 0))/ (a(U » V)) are weaker for the Coulomb interactions since the 
spatial correlations are smaller and do not appear for the short range interactions where no 
correlations exist. 

In order to consider the relevance of the above results to an experimental situation 
one must determine the ration U/V. For quantum dots U/V > 3 is realized, and strong 
evidence of the intermediate range order appears in the behavior of the ground state energy 
fluctuations ||13|| . Therefore, one might expect that these systems will show a strong decrease 
in their polarizability beyond the classical value due to the appearance of intermediate range 
order. 
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FIG. 1. The ratio (a(U = 0))/(a(U)) for short range interactions. The numerical results for 
N = 8 are represented by circles and for N = 4 by squares. The weak interaction regime is plotted 
in (a), while in (b) a larger range of interactions is presented. 
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FIG. 3. The two point correlation C(r = Vbs) for (a) short range interactions, (b) long range 
interactions. Circles correspond to N = 8 and squares to N = 4. 
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